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1.  Technical  Summary 

1.1  Technical  Problem 

The  objective  of  the  research  of  this  contract  is  to  find  temperature 
compensated  materials  for  use  in  surface  acoustic  wave  (SAW)  signal  processing 
devices,  i.e.  materials  with  large  electromechanical  coupling,  low  ultrasonic 
attenuation  and  a vanishingly  small  temperature  coefficient  of  the  delay  time. 

The  electromechanical  coupling  factor  should  be  substantially  larger  than  for 
cx-quartz,  which  is  presently  used  in  temperature  compensated  SAW  devices. 

1.2  Methodological  Approach 

The  research  consists  of:  (A)  both  exploratory  and  systematic  crystal 

growth  studies  on  a variety  of  materials  which  are  expected  to  be  temperature 
compensated  for  bulk  waves  and  which  have  been  selected  earlier  under  AFCRL 
Contract  F19628-73-C-108  on  the  basis  of  certain  heuristic  criteria  (Barsch 
and  Newnham,  1975),  and  (B)  measurements  of  the  single  crystal  elastic  and  thermo- 
elastic properties  of  the  above  grown  crystals  to  determine  whether  they  possess 
temperature  compensated  crystallographic  directions  for  bulk  waves,  and 
measurements  of  their  piezoelectric  and  dielectric  constants  and  their  corres- 
ponding temperature  coefficients  to  check  the  suitability  of  these  materials  for 
surface  wave  device  applications. 

1.3  Technical  Results 

1.3.1  Crystal  Growth 

A.  Crystal  Growth:  The  crystal  growth  effort  centered  on  L^SiO^, 

Ba2Si2TiOg,  and  Pb2KNb,-0^.  Both  Bridgman  and  Czochralski  crystal  growth 
experiments  were  performed  on  L^SiOy  but  the  boules  obtained  were 
polyc.rystalline.  The  high  temperature  behavior  of  this  phase  was  examined  by 
DTA  and  high  temperature  x-ray  techniques  and  indicated  no  major  phase  changes 
between  room  temperature  and  the  melting  point  of  L^SiO^.  Various  shielding 


arrangements  and  temperature  stability  of  the  growth  system  were  examined  to 
determine  their  effects  on  the  quality  of  the  boules.  Boules  of  Ba^SijXiOg 
were  pulled  both  from  stoichiometric  melts  and  from  TH^-rich  melts.  However, 
good  quality  crystals  were  not  obtained  in  either  case.  The  cracking  problem 
in  pulled  Pb2KNo5015  crystals  persisted,  but  was  reduced  by  varying  the  lead 
content  of  the  melt  and  decreasing  the  diameter  of  the  boules.  Several  single 
crystal  pieces  barely  large  enough  for  measurements  were  obtained. 

1.3.2  Measurement  of  Elastic,  Thermoelastic,  Piezoelectric 
and  Dielectric  Properties 

The  theoretical  equations  required  for  the  experimental  determination  of 
the  complete  set  of  piezoelectric  constants  from  X-ray  measurements  have  been 
derived.  For  this  purpose  the  expression  for  the  quantity  (30/0E),  where  0 
denotes  the  Bragg  angle  and  E the  magnitude  of  an  applied  electric  field  is 
calculated  as  a function  of  the  field  direction  and  the  reflecting  lattice  plane 
normal.  For  all  20  crystal  classes  exhibiting  the  piezoelectric  effect  explicit 
expressions  are  given  for  the  longitudinal  and  transverse  piezoelectric  effect, 
corresponding  to  parallel-field  and  perpendicular-field  reflection,  respectively. 
For  the  19  piezoelectric  classes  of  the  monoclinic,  orthorhombic,  tetragonal, 
trigonal,  hexagonal  and  cubic  systems  explicit  expressions  for  (90/0E)  in  terms 
of  the  Miller  indices  of  the  reflecting  planes  ar>  given  for  the  simplest  crystal 
cuts  with  respect  to  the  symmetry  elements  present. 

The  on-diagonal  elastic  constants  c^,  c^,  c^  and  and  their 

temperature  coefficients  have  been  measured  ultrasonically  for  lead  potassium 
niobate,  Pb^b^O^.  The  elastic  constant  values  differ  considerably  from  results 
given  by  Yamada.  The  temperature  coefficient  of  the  longitudinal  modulus  c^  is 
found  to  be  positive,  substantiating  earlier  indications  that  lead  potassium 
niobate  should  have  temperature  compensated  cuts  and  therefore  could  be  a 
superior  substitute  for  a-quartz  in  SAW  devices. 
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1.4  DoD  Implicati  t 

For  one  of  the  materials  investigated,  lead  potassium  niobate,  one  may 
expect  the  existence  of  temperature  compensated  cuts  for  bulk  and  surface  waves, 
with  substantially  larger  electromechanical  coupling  than  for  a-quartz.  Thus 
by  replacing  quartz  as  a substrate  material  in  surface  acoustic  wave  (SAW)  devices 
by  lead  potassium  niobate  insertion  losses  can  be  reduced  and  the  operating 
frequency  and/or  bandwidth  can  be  increased.  In  this  manner  the  efficiency, 
reliability  and  capability  of  military  communications  and  Radar  systems  utilizing 
SAW  signal  processing  devices,  such  as  multichannel  communications,  secure  anti- 
jam communications  for  satellites,  miniature  aivonics  and  electromagnetic  counter 
measures,  can  be  improved. 

1.5  Implications  and  Further  Research 

It  has  been  demonstrated  that  the  search  for  new  temperature  compensated 
materials  with  properties  superior  to  those  of  a-quartz  through  the  approach  used 
under  the  present  contract  can  be  successful.  One  may  therefore  hope  that  a 
continued  systematic  search  for  new  temperature  compensated  materials  under  the 
present  contract  may,  even  with  the  low  funding  level,  eventually  lead  to 
the  discovery  of  additional,  perhaps  even  more  suitable,  materials.  To  this  end 
continued  crystal  growth  efforts  are  required  to  obtain  suitable  single  crystals 
for  the  physical  property  measurements,  which  are  necessary  to  assess  the  use  of 
a given  material  for  SAW  device  applications. 

1.6  Special  Comments 


No  special  comments  are  offered  at  this  time. 


-4- 


•Sf' 


2.  Crystal  Growth  Results 
2,1  Lithium  Silicate,  Li^SiO^ 

Polycrystalline  starting  material  of  this  phase  was  synthesized  at 
900°C  from  Li^CO^  and  SiC^.  Attempts  were  made  to  grow  the  lithium  metasilicate 
from  the  melt  by  both  Bridgman  and  Czo chralski  techniques,  with  the  latter 
appearing  most  promising.  Milky,  polycrystalline  boules  were  all  that  could  be 
obtained  from  the  Bridgman  directional  solidification  experiments,  even  when 
slow  growth  rates  of  0.8  mm,  hr  were  used.  The  temperature  gradient,  at  the  melt- 
crystal  interface  was  about  25  C°/inch  in  these  experiments. 

Czochralski  crystal  pulling  experiments  on  L^SiO^  were  performed  with 
two  different  pieces  of  equipment.  The  first  ones  were  with  an  Arthur  D.  Little 
model  MP  puller  heated  by  an  rf  induction  generator.  Since  cracking  of  the 
boules  persisted  even  after  using  various  shielding  and  insulation  schemes,  we 
decided  to  use  an  NRC  puller  equipped  with  a very  stable  resistance  furnace 
built  for  this  puller.  The  central  furnace  tube  is  wound  with  Pt-10%  Rh,  and 
this  is  surrounded  by  a larger  tube  wound  with  Kanthal.  The  furnace  is  equipped 
with  a proportional  controller  and  a Data-Trak  for  programmable  heating  and 
cooling  cycles.  A special  "cover"  for  shielding  above  the  crucible  was  fabri- 
cated from  an  inverted  kyanite  crucible  fitted  with  a side  arm  and  window  for 
viewing.  The  temperature  gradients  for  the  first  three  inches  above  the  melt 
varied  with  shielding  ciesigns  from  70° /inch  to  200° /inch.  This  cover,  the 
stability  of  the  furnace  (less  than  ± 1/2°  fluctuations),  and  slow  pulling  rates 
of  about  2 ram/hr  greatly  increased  the  perfection  of  the  boules,  but  they 
remained  polycrystalline.  Nucleating  a crystal  on  the  platinum  extension  to  the 
pulling  rod  and  the  necking  of  crystals  has  been  a problem  because  of  (i)  viewing 
difficulties,  (ii)  the  good  shielding  and  resultant  reduced  temperature  gradients 
above  the  melt  surface,  and  (iii)  the  extreme  sensitivity  of  the  crystal  diameter 
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to  the  furnace  temperature.  These  problems  should  be  eliminated  with  the  use 
of  seeds  cut  from  the  polycrystalline  boules  obtained,  and  a newly  designed 
viewing  port. 

The  possibility  of  a phase  transition  causing  the  cracking  in  Li^SiO^ 
was  examined  by  DTA  and  high  temperature  x-ray  experiments.  The  DTA  experiments 
showed  a small  endothermic  peak  at  about  1030°  and  then  the  melting  point  at 
1200°C.  High  temperature  x-ray  diffraction  patterns  taken  at  835°,  1060°,  and 
l\Sba  showed  only  the  expected  thermal  expansion  of  the  orthorhombic  unit  cell; 
no  new  peaks  were  observed.  These  high  temperature  patterns  are  now  being 
examined  to  determine  if  anisotropy  in  thermal  expansion  could  be  great  enough 
to  contribute  significantly  to  the  cracking  problem. 

2.2  Fresnoite,  Ba^Si^TiOg 

In  previous  growth  experiments  on  I^S^TiOg,  its  high  melting  temperature 
of  1360°C  was  difficult  to  maintain  in  the  center  of  the  insulated  platinum 
crucible  without  melting  the  crucible.  This  problem  was  overcome  by  various 
heat  shielding  arrangements  around  the  crucible  as  well  as  reflectors  and  shields 
above  the  crucible,  but  only  poor  quality,  cracked  crystal  boules  were  obtained. 
To  circumvent  the  problem  in  another  way,  flux-pulling  experiments  were  attempted 
since  the  temperature  requirements  are  not  so  severe.  Very  small  crystals  of 
fresnoite  have  been  grown  by  slow  cooling  of  Ti02~rich  melts  at  the  National 
Bureau  of  Standards,  and  our  own  DTA  experiments  indicated  that  additions  of 
TiC>2  to  stoichiometric  fresnoite  lowers  the  melting  temperature  significantly. 
Therefore,  we  carried  out  growth  experiments  after  adding  various  amounts  of 
Ti02  ranging  from  5 to  30  wt.%  to  fresnoite.  In  all  instances  Ba2Si2TiOg  could 
be  pulled  from  these  molten  solutions,  but  all  the  crystal  boules  contained 
numerous  cracks  and  appeared  more  and  more  milky  with  increasing  amounts  of  Ti02 
in  the  melt.  With  30  wt.%  excess  Ti02  the  melting  point  was  lowered  about  150° 
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so  that  a boule  could  be  grown  without  extraneous  heat  shields,  but  the  boules 
were  still  cracked.  A thin  slab  of  crystal  observed  under  the  polarizing 
microscope  showed  appreciable  amounts  of  Ti02  flux  intergrown  with  the  fresnoite. 

Even  with  slow  pulling  rates  of  about  2 ram/hr,  the  TiO^  flux  inclusions  are  still 
a serious  problem. 

2,3  Lead  Potassium  Niobate,  Pb^KNb^O^ 

Several  final  Czochralski  growth  experiments  were  carried  out  on  Pb2KNb5015. 

By  varying  the  content  of  lead  oxide  in  the  melt,  smull  boules  of  approximately 
10  to  15  mm  in  diaiuoter  were  grown  in  which  the  cracking  was  reduced.  Several 
oriented  and  parallel-epiped-shaped  specimens  were  prepared  for  the  ultrasonic 
measurements.  When  viewed  under  a polarization  microscope  one  specimen  appears 
to  be  homogeneous  and  free  of  internal  strains  and  is  therefore  likely  to  consist 
of  a single  domain. 

3,  X-Ray  Determination  of  Piezoelectric  Constants 
3.1  Introduction 

Fo*  the  calculation  of  the  electromechanical  coupling  factor  the  complete 
set  of  piezoelectric  constants  has  to  be  known.  In  the  present  investigation  the 
piezoelectric  constants  are  being  determined  by  measuring  by  means  of  x-rays  the 
elastic  strain  which  is  induced  in  a crystal  by  an  applied  electric  field.  This 
method  was  first  used  by  Bhalla  Bose,  White  and  Cross  (1971)  to  measure  the 
piezoelectric  constant  d^  of  a-quartz.  These  authors  found  good  agreement  with 
earlier  data  obtained  with  other  methods.  As  pointed  out  by  Bhalia,  Bose,  White 
and  Cross  (1971),  the  x-ray  method  has  several  advantages  over  other  methods.  Among 
these  are  its  ready  applicability  to  small  crystals  (about  lramx  1 mm  x 0.1  mm)  and 
the  possibility  to  distinguish  spontaneous  and  induced  3trains  in  ferroelectric 
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crystals  and  to  determine  the  piezoelectric  constants  of  individual  domains. 
In  addition,  the  x-ray  method  avoids  r,ome  of  the  potential  difficulties 
associated  with  other  methods,  such  s s the  occurrence  of  mode  coupling  in 
the  widely  used  resonance-antiresonance  method  (see,  e.g.,  Mason,  1950). 

The  theoretical  equations  for  the  x-ray  dete-rmination  of  the  complete 
set  of  the  piezoelectric  constants  for  crystals  of  all  twenty  crystal  classes 
which  exhibit  the  piezoelectric  effect  have  been  derived.  In  the  following 
the  derivation  and  the  results  will  be.  presented,  and  a numerical  example 
will  be  discussed. 

3.2  Theoretical  Derivation 

Bragg's  form  of  the  considition  for  constructive  reflection  of  an 
incident  x-ray  beam  by  a lattice  plane  with  Miller  indices  hk£  and  inter- 
planar  spacing  d. , ^ J s given  by 

nA 


sin  0 = 


2d 


(1) 


hkl 


where 


= G 


hk£ 


hkA 


1/2 


(2) 


and 


t,  = ha*  + kb*  + He* 
hki 


(3) 


denotes  a reciprocal  lattice  vector,  a*,  b*,  c*  are  the  base  vectors  of 
the  reciprocal  lattice  defined  by  a • a*  ■ 1,  a • b*  ■ 0,  etc.,  where  a,  b, 
c are  the  unit  cell  vectors  of  the  direct  space  lattice.  If,  as  in  equation 


(1),  the  order  of  the  reflection  appears  explicitly  the  Miller  indices  must 
be  assumed  to  consist  of  coprime  integers  (i.e.  of  integers  without  a common 
factor  except  unity). 


By  applying  an  electric  field  E with  components  ■ Eot^,  (c*r  ar e the 
direction  cosines,  and  r = 1,2,3  denotes  three  cartesian  cr  jrdinate  axes) 
to  a crystal,  an  elastic  strain  = 1,2,3)  is  induced  through  the 

converse  piezoelectric  effect  according  to  (Nye,  1957) 


E 


r 


(A) 


d . denotes  the  third  rank  tensor  of  the  piezoelectric  strain  constants, 

rij 

and  the  summation  convention  is  used  here  and  subsequently  (i.e.  summation 

over  the  three  values  1,2,3  is  implied  for  every  pair  of  identical  indices). 

Of  course,  the  components  Er  of  the  electric  field  E and  of  the  strain 

tensor  ^ , and  the  set  of  piezoelectric  constants  drij  must  be  referred  to 

a common  cartesian  coordinate  system  with  axes,  x^,  x^,  x^. 

In  order  to  obtain  the  dependence  of  the  reciprocal  interplanar  distance 

d b on  strain  E,  consider  the  matrices  A and  A*  composed  of  the  unit  cell 
hkx,  ij  fS  35 

vectors  of  the  direct  and  reciprocal  lattices,  respectively,  according  to 


' ax  bx  ' 

[ ci 

b* 

C1  ' 

a2  b2  C2 

, A*  = 

a* 

2 

b*2 

2 

a3  b3  c3 

.‘*3 

C* 

c3 

(5) 


The  relations  a • a*  - 1,  a • b*  = 0 etc.  can  then  be  written  in  matrix 


notation  as 


I 
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£ A*T  = aV  - 1 (6) 

• v <*»■  'V  A> 

where  the  superscipt  T denotes  the  transpose  matrix,  and  1 the  unit 
matrix. 

In  the  approximation  of  the  linear  theory  of  elasticity  under  the 
influence  of  a strain  e the  base  vectors  change  according  to  a(e)  * (1  + e)a(O), 

~ » ~ Tl 

etc.,  where  e denotes  the  strain  tensor  e . . in  matrix  notation.  Therefore, 

~ ij  ’ 

the  direct  and  reciprocal  base  vector  matrices  change  according  to 


A(e)  = 

IS  ^ 

(1  + e)A(O) 

(7a) 

A*(e)  = 

(1  + c)-1A*(0)  - (1  - e)A*(0) 

.CXI  JvC  ***  rV 

(7b) 

In  the  spirit  of  the  linear  theory  of  elasticity,  higher  than  first  powers 
of  j:  have  been  neglected  in  the  last  step  of  equ.  (7b). 

According  to  (7b)  a reciprocal  lattice  vector  changes  according 

to 


Therefore, 


(8) 


and 

Ghk;(e>  ■ Ghw(0)(1  - <s  • 


(10) 


where 
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N - 


’hk  l 


’lik£, 


(ID 


denotes  the  normal  of  the  lattice  plane  (hk£) . 

In  view  of  equsi.  (1),  (2),  (4)  and  (10)  the  diffraction  angle  0 depends  on  the 
magnitude  E of  the  electric  field,  0 a 0(e(E)),  so  that  one  obtains  by 
differentiation  (summation  convention!) 


[din  sin  O') 

(Mn  Ghk*i 

rd£.  .•> 

M 

r3E 

r 

l 3E  J 

l J 

3E 

r 

3E 

or 


cot  0 


30 

3E 


a 


Vj 


(12) 


This  relation  is  the  generalized  form  of  equ.  (6)  of  the  paper  by  Bhalla, 

Bose,  White  and  Cross  (1971).  For  the  special  case  that  both  the  direction 

of  the  electric  field  and  the  reflecting  lattice  plane  normal  are  in  the 

x1  direction,  - 6lr,  N±  - 6^,  N » 5^  (6^  denotes  the  Kronecker  symbol, 

i.e.  6 *1  for  m ■ n = 1,2,3,  and  6 ■ 0 for  m j4  n) , the  above  equ.  (12) 

mn  mn 

becomes  equivalent  to  equ.  (6)  of  these  authors. 

Equ.  (12)  provides  the  theoretical  basis  for  the  experimental  determi- 
nation of  the  complete  set  of  piezoelectric  constants  of  any  piezoelectric 
material.  In  order  to  determine  the  complete  set  of  piezoelectric  constants 
the  dependence  of  the  reflection  angle  0 on  the  magnitude  of  the  electric 

•f 

field  must  be  measured  for  a set  of  combinations  of  directions  a,N  equal  in 
number  to,  or  if  redundancy  checks  are  desired,  larger  than  the  number  of 
independent  piezoelectric  constants  corresponding  to  the  crystal  class  of 
the  material  under  investigation. 

In  order  to  eliminate  electrostriction  the  measurements  should  be  carried 
out  with  field  reversal  which  changes  the  sign  of  equ.  (12),  but  not  of  the 
quadratic  electrostriction  effect. 
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Application  of  the  electric  field  requires  a set  of  thin  platelets 
of  the  piezoelectric  materials  with  plane  parallel  electroded  faces,  so  that 
the  direction  of  one  of  the  two  face  normals  of  opposite  directions  coincides 
with  the  direction  of  3 - (a^c^)  of  the  electric  field.  In  practice,  the 
direction  of  the  reflecting  lattice  planes  is  restricted  to  the  two  cases 
N | j a and  N ]_  a.  While  the  former  case,  corresponding  to  the  reflecting 
lattice  planes  parallel  to  the  platelet  faces,  utilizes  the  longitudinal 
piezoelectric  effect  and  is  more  convenient  experimentally,  it  will  be  shown 
that,  with  the  exception  of  four  crystal  classes  (6,  6m2,  23,  43m)  the  complete 
set  of  piezoelectric  constants  cannot  be  determined  in  this  manner.  On  the 
other  hand,  the  second  configuration,  which  utilizes  the  transverse  piezo- 
electric effect  and  corresponds  to  the  electric  field  lying  in  the  reflecting 
planes,  requires  larger  experimental  effort  to  eliminate  distortion  of  the 
electric  field  near  the  side  faces  of  the  sample,  but  permits  to  determine  the 
complete  set  of  piezoelectric  constants  for  all  twenty  piezoelectric  crystal 

classes. 

In  the  following  section  explicit  expressions  of  the  right-hand-  ide 
of  equ.  (12)  will  be  given  for  all  twenty  piezoelectric  crystal  classes  and 

for  both  experimental  configurations. 

3.3  Application  to  Piezoelectric  Crystal  Classes 

Denoting  for  convenience  the  RHS  of  equ.  (12)  by  L and  T for  parallel 
field  and  perpendicular-field  reflection,  corresponding  to  the  longitudinal 
and  transverse  piezoelectric  effects,  respectively,  equ.  (12)  may  be  rewritten 


as 
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- cot  0 


fas)  JL  <“  M s> 

,3e J " \ T (a  J.  it) 


where 

L"*Wjdrij 
T - d (a  • ft  - 0) 


(13) 


(14a) 

(14b) 


For  the  numerical  application  of  these  equations  the  form  of  the 
piezoelectric  constant  tensor  corresponding  to  a particular  crystal  class 
must  be  used,  and  both  the  lattice  plaae  normal  N and  the  direction  a of  the 
electric  field  must  be  expressed  in  terms  of  the  Miller  indices.  For  the 
components  of  the  unit  vector  N the  corresponding  relation  is  obtained  from 
equs.  (3)  and  (11)  and  is  for  the  general  triclinic  case  given  by  the 
familiar  expression  (International  Tables  for  X-Ray  Crystallography  [1959]) 
(i  - 1,2,3): 


ha*+kb*+)tc* 


,2  *2  2 *2  -2  *2  .**  * „*  * * **  * 1/2 
[h  a +k  b +SL  c + 2k£b  c cos  a +2h£.a  c cos  3 +2hka  b cos  y ] 


(15) 


* * * , * * * 

Here  a , b , c denotes  the  lattice  constants,  and  a , 3 , Y the  angles 

-+■*  ■+■*  -+■*  * * 
between  the  unit  cell  vectors  a , b , c of  the  reciprocal  lattice,  a^,  b^, 

* -+•*  -+•*  -+* 
c^  (i=l,2,3)  denotes  the  cartesian  coordinates  of  a , b , c , respectively, 

in  the  same  coordinate  system  to  which  the  piezi  electric  constant  tensor  is 

referred. 

For  parallel  field  reflection  the  unit  vector  a is  determined  by  the 

->■  -f  -► 

condition  a ■ N.  For  perpendicular  field  reflection  a is  constrained  by  the 

->•  -f- 

condition  a*N  * 0.  Thus  for  a given  lattice  plane  (hkH)  the  vector  a has  one 
degree  of  freedom,  which  may  be  conveniently  described  by  the  angle  <f>  according 
to 
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j£  ■ cos  <}>  U + sin  $ ^ 


(16) 


where  U is  a unit  vector  perpendicular  to  N lying  in  the  plane  generated  by  the 
unit  vector  ^ in  the  x3  direction,  and  by  the  lattice  plane  normal  N,  and  V 
is  perpendicular  to  both  U and  N.  The  vectors  U and  V are  then  given  by 


U = (-(e3*N)N  + e3)/(N3  + N2)1/2 


V = U x N 


(17a) 


(17b) 


Inserting  (17a)  and  (17b)  into  (16),  the  components  a.  of  the  unit  vector  a 
may  be  expressed  in  terms  of  the  components  of  the  lattice  plane  normal 

and  of  the  angle  <j>  according  to 

* - (N^N3  cos  4>  + N2  sin  <£)/  (N^  + N^)  ^ ^ 

a2  = - (N2N3  cos  4)  - Nx  sin  <fr)/(N*  + N2)l/2  (18b) 

2 2. 1/2  * (18c) 

a3  = ^N1  + ^2  C°S  ^ 

As  will  be  discussed  below,  lor  most  crystal  classes  convenient  choices 
of  the  angle  * are  possible  ll  at  least  one  of  the  Miller  indices  is  tero. 

In  these  cases  the  vector  a lies  along  symmetry  directions  of  the  crystal, 
and  a considerable  simplillcation  of  the  equations  given  below  results. 

By  using  the  form  of  the  piezoelectric  tensor  corresponding  to  the 
individual  crystal  classes  as  given,  for  example,  by  Xye  (1957)  the  expressions 
L and  T defined  in  (Ha)  and  (14b)  may  be  explicitly  written  out  as  given  below 
for  the  twenty  piezoelectric  crystal  classes.  In  the  following,  the  piezo 
electric  constants  will  be  expressed  in  Voigt  notation,  that  is,  the  index 
pair  ij  and  drlJ  (with  r.i.J  ■ 1.2.3)  is  replaced  by  a single  index 
p - 1, 2,3:4; 5(6  corresponding  to  ij  = 11,22,33)23  and  32  il3  and  31)12  and 
21,  respectively,  and  a factor  of  two  is  introduced  for  p - 4,5,6,  such  that 


-14- 


d = dfij  for  p ■ 1,2,3,  and  df^  ■ 2d^j  ^or  ^ “ 4,5,6.  The  following 
abbreviations  are  introduced  in  order  to  describe  the  coefficients  of  the 
piezoelectric  constants  occurring  in  the  expression  for  T for  the  tetragonal, 
trigonal,  hexagonal  and  cubic  systems: 


22  22  1/2 
A+  *=  [-2N1N2N3  cos  p + (N1-N2)sin  $]N3/ (N^+Np 


2 2 1/2 

A_  = - N3(N^+N2)  ' sin  <f> 

222  1/2 

B+  + _ N3(N^+N2)  cos  p 

22  22  1/2 
B_  « - [(N^-N2)N3  cos  4>  + 2N1N2  sin  ^^/(Np-Np 

C+  - (N^+N2)3/2  cos  p 

c_  - (n2-n2)(n2+n2)1/2  cos  ^ 

22  22  221/2 

Dx  - - [(N^-3N2)N1N3  cos  p + (3N1-N2)N2  sin  (jiJ/O^+Np 

22  22  221/2 
D2  - [(3N^-N2)N2N3  cos  p - (N^Npt^  sin  <J)]  / (N1+N2> 

22  22  22  1/2 

E = [(3(N3+N2)-2)N1N2  cos  p + (N1-N2)N;j  sin  ($> ] / (N1+N2) 


222  1/2 

F = N3(Nx+N2)  ' cos  p 


2 2 1/2 

G «=  N1N2(N1+N2)  cos  <j> 


(19a) 

(19b) 

(19c) 

(19d) 

(19e.) 
(19f ) 

(19g) 

(19h) 

(19i) 

(19  j ) 

(19k) 


The  description  of  the  direction  a of  the  electric  field  in  terms  of 
the  angle  (}>  as  given  in  equs.  (18),  and  the  U9e  of  the  ^-dependent  quantities 

0 -f 

defined  in  equs.  (19a)  to  (19k)  are  limited  to  the  case  N f e3>  For  N = e3 
the  angle  p is  not  defined.  However,  the  special  case  N3  = e3  is  of  interest 
only  for  the  triclinic  and  monoclinic  systems,  and  the  quantity  T defined  in 
equ.  (14b)  becomes  in  this  case 
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T = ald13  + a2d23 


(20) 


On  the  other  hand,  for  the  tetragonal,  trigonal,  hexagonal  and  cubic  systems 
the  quantity  T is  identically  zero  for  N = e^,  so  that  this  case  need  not 
be  considered  explicitly. 

1.  Trlclinic  System 
Class  l(Cj.) 

(18  independent  constants:  d^,  d12»  di;j,  d.^,  d^,  d1&,  d21>  d22>  d23,  d^, 

d25*  d26’  d31 ' d32 * d33’  d34’  d35’  d36) 


"JdU  + N1 

N2(d12  d 

■ d26>  + V 

'3<d13  + 

d35) 

NlN2(d21  + 

d16}  + 

»2d22  + >¥ 

l3(d23  + 

d34> 

NlN3(d31  + 

d15)  + 

N2N3(d32  + 

d2<>  + 

N3d33 

NlN2N3(d14 

+ d25  ‘ 

' d36) 

°lNldll  + 

alN2d12 

+ 0lN23d13 

“2NLd21  + 

a2N2d22 

+ C-2N3d23 

“3Nld31  + 

„2, 

a3N2d32 

+ a3N3d33 

aiN2N3d14 

+ a^N^N 

3d15  + alNlN2d16 

a2N2N3d24 

+ c^NjN 

3d25  + a2N  l1 

V26 

a3N2N3d34 

+ cijNjN 

3d  35  + a3Nl' 

^2d36 

(21a) 


( 21b) 


2.  Monoclinic  System 
Class  2(C 2) 

(8  independent  constants:  d^,  ^16*  d21’  d22’  d23’  d25’  d34*  d36^ 

L - N^d22  + »iN2(d21  + d16)  + + d34) 


+ »l»3W3<a»  + d25  + d36) 


(22a) 
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T = a2Nld21  + a2h2d22  + a2N3d23 


+ a1N2N3d14  + aiNlN2d16  + a2NlN3d25 


+ a3N2N3d34  + a3N1N2d36 


(22b) 


Class  m(C  ) 


(10  independent  constants:  d^,  d^,  d^,  d^,  d^>  d^,  d^,  d32>  d.^,  d33> 


L - Nddu  + NlN2(d12  + d26)  + NlN3(d13  + d35) 
+ NjN3(d31  + d15)  + N22N3(d32  + d24)  + N^3 

T ‘ alNldll  + aiN2d12  + alN3d13 
+ a3Nxd31  + a3N2d32  + a3N3d33 


(23a) 


+ a1N1N3d15  + a2N2N3d24  + a2NlN2d26  + a3NlN3d35 


(23b) 


3.  Orthorhombic  System 
Class  222 (D  ) 

(3  independent  constants:  d^,  d23>  d^ 

L - + d25  ♦ d36) 


T = ot1N2N3di^  + a2NiN3d25  + a3NlN2d36 


(24a) 

(24b) 


Class  mm2(C2v) 

(5  independent  constants:  d^,  d^,  d^,  d32»  d33) 


L = N?N0(d,,  + d1K)  + N2N3(d32  + d9/)  + N2d 


13  31  15 


24' 


3 33 


T * »3N1d31  + a3N2d32  + a3N3d33 


(25a) 


+ a.N.N-d.,.  + a0N0N,d„, 
1 1 3 15  2 2 3 24 


(25b) 
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4.  Tetragonal  System 
Class  4(C,) 

(4  Independent  constants:  d^,  d^,  d^,  d33) 

1.  - (N‘  + t*2>N3<d31  + d15>  + N^d33 

(26a) 

1 * C+d31  + Fd33  + A-d14  + B+d15 

(26b) 

Class  4(S, ) 

(4  independent  con., rants:  d^  d^,  d^,  d^) 

L ' <N1  - N2,N3(d31  + “is’  + W3<2dU  + d36> 

(27a) 

I - C_d31  + A,du  + B..dls  + Gd36 

(27b) 

Class  422(b,) 

(1  independent  constant:  d,^) 

L = 0 

(28a) 

T " A-d14 

(28b) 

Class  4nnn(C,..) 

(3  independent  constants:  d^,  d33) 

L - ("J  + N2>N3(d31  + d15>  + K3d33 

(29a) 

1 ■ C+d31  + Fd33  + B+d15 

(29b) 

Class  A2m(D^) 

(2  Independent,  constants:  d.^,  d^) 


L = N1N2N3(2du  + d3(.)  (30a) 

36 


T = A,d..  + Gd 
+ 14 


(30b) 
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5.  TrlRonal  System 
Class  3(C^) 

(6  independent  constants: 


L - (Nj  - 3N2)N1dn 


+ (Nj  + N^)N3(d31 


T = Dldll  + °2d22  + 


dll*  d14 ’ d15’  d22 * d31*  d33) 
(3Ni  " N2'N2d22  + N3d33 

+ *15* 

C+d31  + FD33  4-  A_du  4 B+d15 


Class  32 (D  ) 

(2  independent  constants:  d^, 

L - (nJ  - 3N2)N1d11 


Dldll  + A-d14 


Class 


_Mc3vI 


(4  independent  constants:  d.^,.,  d^,  ri31>  d33) 

L ■ - (3NJ  - N^)N2dss  + <»J  + N*)N3(d31  + djj)  + Ndd33 


V22  + C+d31  + FdJ3  + S+d15 


6.  Hexagonal  System 
Class  6(C^) 

(4  independent  constants:  d^,  d^,  d^,  d33> 

L ' <N1  + N2>N3<d31  + d15>  + N3d33 
T " c+d31  + Fd33  + A-d14  + Ml5 


Class 


(2  independent  constants:  d^,  d ) 

L - (Nj  - - (3N^  - 


T ” °ldll  + D2d22 


(31a) 

(31b) 


(32a) 

(32b) 


(33a) 

(33b) 


(34a) 

(34b) 

(35a) 

(35b) 


! 
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Class  622  (D^.) 

(1  independent  constant:  d^) 


Ii  = 0 

(36a) 

T - A-du 

Class  6mm (C^) 

(3  independent  constants:  d^,.«  u3i*  ^33^ 

(36b) 

L - <Nj  + N2)N3(d31  + dl5)  + N3d33 

(37a) 

T " C+d31  + Fd33  + Mis 

Class  6m2(D-j^) 

(1  independent  constant:  d.^) 

(37b) 

L - (Nj  - 3N^)N1d11 

(38a) 

T - Dldn 

Classes  23 (T)  and  43m(Tj) 

(1  independent  constant:  d^) 

(38b) 

L - 3N1N2N3d14 

(39a) 

T = Ed14 

(39b) 

As  mentioned  above,  equations  (26b),  (27b),  etc.,  to  (39b)  are  valid  only 
for  N i ty  and  for  N = the  quantity  T in  these  equations  is  identically 
zero. 

Equation.  (21a)  to  (39a)  represent  the  most  general  form  of  the 
expressions  for  the  longitudinal  piezoelectric  effect  for  the  20  piezoelectric 
crystal  classes.  Each  lattice  plane  (hk£)  is  characterized  by  its  effective 
piezoelectric  coefficient  L = L^,  which  is  for  triclinic  symmetry  obtained 
by  inserting  cqu.  (15)  into  cqu.  (21a),  and  for  the  crystal  classes  of  higher 
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symmetry  by  inserting  the  simplified  form  of  equ.  (15)  appropriate  for 
the  particular  crystal  symmetry  into  the  corresponding  equs.  (22a)  to  (39a). 

For  each  lattice  plane  (hkJl)  the  electric  field  coefficient  -cot  0 ( 3 C / 3 K ) 
of  the  Bragg  angle  0 = Oj ^ gives  according  to  equ.  (13)  an  experimental 
value  for  the  corresponding  linear  combination  of  piezoelectric  constants. 

While  no  further  simplification  is  possible  for  triclinic  symmetry,  it  will 
be  shown  below  that  for  the  crystal  classes  of  higher  symmetry  the  expressions 
for  L given  by  equs.  (22a)  to  (39a)  may  become  much  simpler  if  one  or  two 
of  the  Miller  indices  are  zero.  However,  in  several  instances  the  coefficients 
of  the  piezoelectric  constants  in  the  expressions  for  L vanish,  unless  all 
three  Miller  indices  are  nonzero.  Thus  for  the  experimental  determination 
of  the  complete  set  of  piezoelectric  constants  the  most  general  form  of 
the  expressions  for  L as  given  in  equs.  (21a)  to  (39a)  has  to  be  used. 

Similarly,  equ-tiens  (21b)  to  (39b)  represent  the  most  general  form 
of  the  expressions  for  the  transverse  piezoelectric  effect  in  terms  of  the 
components  of  the  reflecting  lattice  plane  normal  N,  the  direction  a of  the 
electric  field,  and  the  coefficients  A+,  A_  etc.  defined  in  equs.  (19a)  to 
(19k).  According  to  equs.  (18a)  to  (18c)  the  components  of  a,  and  according 
to  equs.  (19a)  to  (19k)  the  coefficients  A+,  A , etc.  depend  on  the  components 
of  N and  on  the  free  parameter  4>  defined  in  equ.  (16).  The  components  of  N, 
in  turn,  depend  on  the  Miller  indices  according  to  equ.  (15).  Thus  all 
coefficients  of  the  piezoelectric  constants  in  the  expressions  for  T depend 
on  the  Miller  indices  in  a rather  involved  manner.  While  for  triclinic 
symmetry  the  equations  have  to  be  used  in  this  form,  a considerable  simplifi- 
cation occurs  if  one  or  two  of  the  Miller  indices  are  zero.  However,  for  the 
two  monoclinic  crystal  classes  2 and  m the  coefficients  of  some  piezoelectric 
constants  (those  of  and  d^  for  class  2,  and  of  for  class  m)  vanish 


-21- 


in  these  special  cases.  Therefore  the  complete  set  of  piezoelectric 
constants  cannot  be  determined  only  from  measurements  of  the  transverse 
piezoelectric  effect  perpendicular  to  these  special  lattice  planes.  In 
these  cases  one  has  to  use  the  general  form  of  equations  (22b)  and  (23b) 
for  the  transverse  piezoelectric  effect  or,  alternatively,  supplement  the 
equations  for  the  transverse  piezoelectric  effect  for  the  special  set  of 
lattice  planes  with  one  Miller  index  zero  by  the  general  equations  (22a)  and 
(23a)  for  the  longitudinal  piezoelectric  effect. 

In  Tables  1 to  5 the  simplified  relations  obtained  from  equs.  (22a) 


and  (22b)  to  (38a)  and  (38t)  are  presented  that  result  for  the  longitudinal 
and  transverse  piezoelectric  coefficients  L and  T in  the  monoclinic,  ortho- 


rhombic, tetragonal,  trigonal  and  hexagonal  classes  if  at  least  one  of  the 
Miller  indices  is  zero.  For  the  monoclinic  and  orthorhombic  classes  the 
general  expressions  for  the  longitudinal  piezoelectric  coefficients  L have 
also  been  included  because  they  are  needed  for  the  experimental  determination 
of  the  complete  set  of  piezoelectric  constants.  For  completeness,  the 
equation  for  L for  the  two  isometric  classes  23  and  43m  is  also 


presented  as  Table  6.  The  explanation  of  Tables  1 to  6 is  given  below. 

The  simpler  equations  presented  in  Tables  1 to  6 are  sufficient  for 
the  experimental  determination  of  the  complete  set  of  piezoelectric  constants 
for  all  crystal  classes,  except  for  triclinic  symmetry.  Therefore  the  general 
form  of  the  equations  (22a)  and  (22b)  to  (39a)  and  (39b)  need  not  be  considered 
in  the  design  of  actual  experiments.  It  may  happen  in  certain  cases,  however, 
that  for  the  lattice  planes  with  at  least  one  Miller  index  zero  that  are 


listed  in  Tables  1 to  6 the  structure  factor  may  be  inconveniently  small, 
or  vanish  exactly.  In  these  cases  and  for  tricllnic  symmetry  one  has  to  resort 


to  the  general  equations  (21)  to  (39) • 


ible  1.  Crystal  plate  normal  a ■ (a.a-a.),  Miller  indices  hkJ. , face  normal  of  reflecting  plane  N “ (N  l^NO  , and 

1 1 J f 

expressions  for  L and  T defined  in  equations  (14a)  and  (14b)  for  monoclinic  system. 

The  ba;e  vectors  b and  c coincide  vi.h  the  Y and  Z axes,  respectively,  and  the  base  vector  a lies  in  the 
X-Z  plane,  forming  an  angle  0 > 90°  with  the  c.  axis. 
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Table  2.  Crystal  plate  normal  a - (a  a^) , Miller  indices  hid,  face  normal  of  reflecting  plane  N - 
and  expressions  for  L and  T defined  in  equations  (14a)  and  (14b)  for  orthorhombic  system. 
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Table  4.  Crystal  plate  normal  a - (a^ctj) , Miller-Bravais  indices  hk£,  for  hexagonal  indexing,  face  normal 
of  reflecting  plane  N = (N^t^)  , and  expressions  L and  T defined  in  equations  (14a)  and  (14b)  for 
trigonal  system. 
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Table  6.  Crystal  plate  normal,  Miller  indices  (hki) , face  normal  of 
reflecting  plane  N = (N^N^),  and  expression  L defined  in 
equation  (14a)  for  isometric  classes  23  and  43m. 
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It  is  apparent  from  equations  (21a),  (22a),  etc.,  to  (39a),  that 
some  of  the  piezoelectric  constants  enter  the  expressions  for  L only  in 
the  form  of  linear  combinations,  so  that  the  complete  set  of  piezoelectric 
constants  cannot  be  determined  from  parallel-field  reflection  measurements 
(a  *»  N)  alone.  Exceptions  are  the  four  classes  6,  6m2,  23,  and  43m,  for 
which  the  quantity  L depends  on  the  individual  piezoelectric  constants  only, 
so  that  they  can  be  determined  from  parallel-field  reflection  alone.  For 
all  other  classes,  perpendicular-field  reflection  measurements  (a^  • N = 0) 
are  required  in  addition  to,  or  instead  of  parallel-field  refJection 
measurements . 

Another  noteworthy  feature  of  the  above  equations  is  that  for  the  six 
classes  4,  422,  3,  32,  6,  622,  the  piezoelectric  constant  d^  does  not  enter 
the  equations  for  L and  therefore  cannot  at  all  be  determined  from  parallel- 
field  reflection  measurements. 

According  to  equations  (21b),  (22b),  etc.,  to  (39b)  the  quantity  T 
depends  in  all  cases  on  all  individual  piezoelectric  constants.  Exceptions 
are  particular  directions  of  N and/or  values  of  <j),  for  which  some  of  the 
coefficients  given  in  equations  (19a)  to  (19k)  may  become  equal,  or  vanish. 
In  general,  however,  it  is  possible  to  determine  the  complete  set  of  piezo- 
electric constants  from  a sufficiently  large  number  of  perpendicular-field 
reflection  measurements  alone.  In  practice,  however,  it  is  desirable  to  use 
as  many  parallel-field  reflection  measurements  as  possible,  because  the 
specimen  preparation  is  simpler  and  experimental  errors  may  be  expected  to 


be  smaller  in  this  case. 
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3.4  Specimen  Orientation 

For  the  measurement  of  the  longitudinal  piezoelectric  effect  crystal 
platelets  with  electroded  pairs  of  (hkJl)  faces  with  normal  1$  ■ a are  required. 
For  the  transverse  piezoelectric  effect  platelets  with  electroded  pairs  of 
faces  with  normal  a perpendicular  to  the  normal  N of  the  reflecting  (hkJl)  face 
must  be  prepared.  In  both  cases  it  is  desirable  to  choose  the  direction  N, 
and  in  the  second  case  also  the  direction  a,  along  high  symmetry  directions 
(parallel  to  rotation  axes,  perpendicular  to  or  within  mirror  planes)  so  as 
to  achieve  accurate  orientation  and  reduce  the  effort  required  for  the  orienta- 
tion and  preparation  of  the  platelets.  However,  since  either  or  both  piezo- 
electric effects  may  vanish  along  some  symmetry  directions  the  requirement 
of  non-vanishing  piezoelectric  effect  eliminates  many  otherwise  suitable 
lattice  planes.  It  is  further  desirable,  that  for  a given  platelet  normal  a 
several  non-equivalent  side  faces  (hk£)  exist,  which  lead  to  different  linear 
combinations  of  piezoelectric  constants  in  the  expression  for  the  transverse 
piezoelectric  coefficient  T.  Ideally,  it  would  be  desirable  to  determine  the 
complete  set  of  piezoelectric  constants  from  one  polygon-shaped  crystal  plate- 
let with  several  side-faces,  so  as  to  require  only  one  pair  of  electroded 
surfaces  for  measuring  the  longitudinal  piezoelectric  effect  in  the  direction 
of  the  plate  normal  a,  and  the  transverse  effect  for  several  non-equivalent 
side  face  normals  N.  Unfortunately,  this  ideal  design  is  possible  only  for 
some  of  the  high  symmetry  crystal  classes.  In  general,  more  than  one  type  of 
crystal  platelet  is  required  for  determining  the  complete  set  of  piezoelectric 
constants.  In  this  case  it  is  desirable  to  perform  the  measurements  mostly, 
or  if  possible,  exclusively  on  the  basis  of  the  longitudinal  effect,  because 
the  corresponding  platelets  require  only  one  pair  of  parallel  faces  and  are 
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therefore  easier  to  prepare  than  platelets  for  the  measurement  of  the 
transverse  effect,  and  because  measurements  of  the  longitudinal  effect  can 
be  done  more  accurately  and  more  easily  than  those  of  the  transverse  effect. 

As  mentioned  above,  this  can  only  be  achieved  for  the  four  classes  6,  6m2, 

23  and  A 3m. 

For  the  triclinic  system  no  symmetry  directions  exist,  so  that  (at 
least)  18  platelets  corresponding  to  18  different  face  normals  a are  required. 
Since  according  to  equation  (21a)  the  longitudinal  piezoelectric  coefficient 
L depends  on  10  linear  combinations  of  piezoelectric  constants,  10  of  the  18 
platelets  may  be  chosen  such  that  their  faces  correspond  to  10  different 
(hki.)  planes.  The  remaining  (minimum  of)  8 platelets  must  be  prepared  for 
the  transverse  piezoelectric  effect  in  such  a manner  that  the  normals  ot  of 
the  pair  of  large  electrodea  surfaces  are  perpendicular  to  8 different  (hk£) 
side  faces.  The  resulting  piezoelectric  coefficients  depend  on  the  orienta- 
tion of  the  crystallographic  axes  with  respect  to  the  cartesian  coordinate 
system  chosen.  It  is  recommended  to  use  the  base  vectors  of  the  direct  and 
reciprocal  lattice  according  to  the  1949  IRE  convention  (Standards  on  Piezo- 
electric Crystals,  1949): 

a = (a  sin  8,  0,  a cos  8) 

b = (-b  sin  a cos  y*,  b sin  a sin  y*,  b cos  a) 

c = (0,  0,  c) 

a*  = (a*  sin  y*,  a*  cos  y*,  0) 

b*  = (0,  b*,  0) 

c*  = (-c*  sin  a*  cos  8.  c*  cos  a*,  c*  sin  a*  sin  8) 


r 
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Here  It  is  (International  Tables  for  X-Ray  Crystallography,  1959): 


a* 


be  sin  a 
V 


etc. , 


V 


2abc  sin 


a+g+Y 

2 


sin 


-q+g+Y 

2 


sin 


q-g+Y 

2 


sin 


q+g-Y 

2 


cos  a*  - (cos  g cos  Y ” cos  a) /sin  g sin  Y » etc. 

It  is  further  q > 90°  , g > 90® 

For  the  symmetry  classes  of  higher  symmetry  the  1949  IRE  convention 
(Standards  on  Piezoelectric  Crystals,  1949)  is  easily  summarized  as  follows. 
The  twofold  axes  of  class  2 (C2) , and  the  mirror  plane  normal  of  class  ra(Cg) 
coincide  with  the  crystallographic  b axis,  which  is  ; aced  along  the  positive 
Y direction  of  a right-handed  cartesian  coordinate  system.  The  c-axis  is 
placed  along  the  positive  Z-direction,  and  the  a-axis  (being  normal  to  the 
b-axis  and  forming  an  angle  g > 90°  with  the  c-axis) is  lying  in  the  XZ-plane. 
For  the  orthorhombic,  tetragonal  and  cubic  systems  the  crystallographic  axes 
coincide  with  the  X,Y  and  Z-axes.  For  the  trigonal  and  hexagonal  systems 
the  threefold  (or  sixfold)  axis  is  parallel  to  the  Z-axis,  and  the  positive 
X-axis  coincides  with  any  of  the  three  secondary  axes  a^  a2>  a3  perpendicular 
to  the  threefold  (or  sixfold)  axes.  The  Y-axis  is  perpendicular  to  the  X- 
and  Z-axes,  so  as  to  form  a right-handed  cartesian  coordinate  system. 


3.5  Discussion  of  Tables  1 to  6 

For  the  remaining  19  piezoelectric  crystal  classes  the  existence 
of  symmetry  elements  reduces  not  only  the  number  of  independent  piezo- 
electric constants,  but  also  leads  to  the  existence  of  high-symmetry 
directions  a with  non-vanishing  longitudinal  and/or  transverse  piezoelectric 
effect,  in  Table  1 to  6 such  pairs  of  high-symmetry  directions  a,  N,  the 
associated  Miller  indices  hkil  and  the  quantitites  L or  T describing  the 


longitudinal  and  transverse  piezoelectric  effect  as  defined  in  equations 
(14a)  and  (14b)  are  listed  for  the  crystal  classes  of  the  monoclinic, 
orthorhombic,  tetragonal,  trigonal,  hexagonal  and  cubic  systems,  respectively. 
In  the  crystal  classes  of  lower  symmetry  the  number  of  high  symmetry 
directions  is  too  small  for  determining  the  complete  set  of  the  piezoelectric 
constants.  In  these  cases  the  expressions  for  both  a and  N and  for  the 
associated  values  of  L and/or  T are  also  given  for  a general  set  of  Miller 


indices . 


The  results  in  Table  1 indicate  that  for  the  class  0^(2)  the  8 

independent  constants  may  be  determined  from  5 crystal  platelets.  From 

*► 

one  platelet  with  orientation  a = (010)  it  is  possible  to  determine  the 
four  constants  d ^ and  d2^»  but  for  the  remaining  four  constants 

^14’  ^16’  d34  3nd  dJ6  °ne  crystal  pdate  each  is  required.  If  the  platelet 
with  orientation  a = (010)  has  3 side  faces  corresponding  to  3 sets  of 
non-equivalent  Miller  indices  h0£  measurement  of  the  longitudinal  effect 
gives  d22»  and  of  the  transverse  effect  perpendicular  to  each  of  the  3 
side  faces  gives  3 linear  combinations  each  of  the  constants  d2^,  d2^  and 
d2y  Of  the  remaining  four  measurements  three  may  be  made  by  using  the 

longitudinal  effect  and  give  (if  the  above  mentioned  constants  d2^,  d22» 

d22  and  d2^  have  been  determined  independently)  the  constants  d^,  ^34*  ^34 

+ d_,.  In  order  to  determine  the  constants  d,  , and  d_.  individually,  one 
36  14  36 

additional  measurement  based  on  the  transverse  effect  is  required. 

For  the  determination  of  the  10  independent  constants  of  class  C^(m) 

9 platelets  with  different  orientations  are  required.  While  the.  platelet 
with  a = (100)  permits  to  determine  d^  through  the  longitudinal  effect, 
and  d^  through  the  transverse  effect,  the  remaining  8 constants  require 


for  their  determination  one  platelet  e ich.  While  5 linear  combinations  of 
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these  constants  (d26,  d13  + d^,  d.^  + d15>  d32  + d^,  d33)  may  be 
determined  from  the  longitudinal  effect,  the  remaining  thre«  independent 
measurements  must  be  based  on  the  transverse  effect. 

As  shown  by  the  data  in  Table  2 the  3 independent  piezoelectric 
constants  of  class  222  have  to  be  determined  from  the  transverse  effect  by 
using  3 different  crystal  platelets  with  faces  normal  to  the  3 coordinate 
axes,  respectively.  The  longitudinal  effect  on  a platelet  with  a pair  of 
general  (hki)  faces  permits  an  independent  measurement  of  the  sum  d17  + 


The  determination  of  the  5 independent  constants  for  class  mm2(C2v) 
is  illustrated  below  with  the  aid  of  a numerical  example. 

The  discussion  of  Tables  3 to  6 for  the  remaining  crystal  classes 
proceeds  along  the  same  lines  as  above  and  need  not  be  explicitly  presented 
here. 


3.6  Numerical  Example 

In  order  to  illustrate  and  facilitate  the  use  of  Tables  1 to  6 the 
expressions  for  L and  T required  for  the  determination  of  the  5 independent 
constants  of  Bi^iNbO^  (crystal  class  mm2)  will  be  given  here.  The  crystal 
structure  of  Bi^iNbO^  has  been  refined  from  x-ray  and  neutron  diffraction 
data  by  Wolfe,  Newnham,  Smith  and  Kay  [1971].  These  authors  list  the 
observed  and  calculated  absolute  values  of  the  structure  factors  of  a large 
number  of  non-equivalent  reflections,  from  which  suitable  high  intensity 
reflections  can  be  selected. 

According  to  Table  2 the  3 constants  d^,  d32  and  d33  can  be  determined 
from  one  platelet  with  a pair  of  electroded  (001)  faces,  if  the  longitudinal 
effect  is  measured  by  Bragg  reflection  on  the  (001)  face,  resulting  in  d(3, 
and  if  the  transverse  effect  is  measured  by  Bragg  reflection  on  two  different 
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wide  faces  of  the  type  (hkO) , which  gives  two  different  linear  combinations 
of  d^  and  d^*  Based  on  the  structure  factors  listed  by  Wolfe,  Newnham, 

Smith  and  Kay  [1971]  we  select  the  three  reflections  0 0 10,  040,  and  400. 

The  corresponding  face  normals  N are  (0,0,1),  (0,1,0)  and  (1,0,0),  respectively. 

o o o 

From  the  lattice  constants  a = 5.431  A,  b = 5.389  A and  c = 25.050  A given 
by  Wolfe,  Newnham,  Smith  and  Kay  [1971]  the  corresponding  Bragg  angles  are 
calculated  for  CuK^  radiation  and  are  listed  in  Table  7,  together  with  the 
values  of  hki,  with  the  structure  factors,  and  with  the  values  of  the  platelet 
normal  a and  the  values  N of  the  three  reflecting  lattice  planes.  In  Figure  la 
the  shape  and  orientation  of  the  platelet  with  faces  corresponding  to  these  three 
normals  and  the  corresponding  Bragg  angles  0^,  and  0^  are  illustrated  schemat- 
ically. Also  listed  in  Table  7 are  the  corresponding  values  of  L and  T obtained 
from  the  last  column  of  Table  2.  According  to  equation  (13)  these  quantities  are 
equal  to  the  electric  field  coefficient  - cot  G^SQ^/SE)  (i  = 1,2,3)  of  the 
associated  Bragg  angle.  Denoting  the  experimental  values  of  the  electric 
field  coefficient  by  M^(i  = 1,2,3)  the  three  piezoelectric  constants  d^, 
d^2  and  d^  are  then  given  by  the  first  three  equalities  listed  in  the  last 
column  of  Table  7.  According  to  Table  2 the  remaining  two  constants  d^  and 
d£4  have  to  be  determined  from  the  transverse  piezoelectric  effect  on  two 
different  platelets  with  side  faces  (h0£)  and  (OkJl)  , with  the  corresponding 
normals  of  the  side  faces  (£0?)  and  (OnC) , respectively,  and  with  normals 
(£0Q  and  (O^n)  of  the  electroded  main  faces.  In  the  fifth  column  of  Table  2 
the  components  (£n(,)  of  these  two  normals  are  expressed  in  terms  of  the  hk&  values. 

On  the  basis  of  the  structure  factors  listed  by  Wolfe,  Newnham,  Smith  and 

Kay  [1971]  we  select  the  reflections  0 2 10  and  2 0 10.  In  Table  7 their  structure 

-V 

factors,  the  normals  N of  the  side  faces  and  the  normals  a of  the  electroded 
main  faces  as  calculated  from  columns  2,  4 and  5 of  Table  2 are  listed 


E 

i 


Miller  indices  hki,  . alculated  structure  factors  |Fj  (in  units  given  by  Wolfe,  Nevnham,  Smith  and  Kay  [19711), 
face  normals  N of  crystal  platelets  representing  reflecting  lattice  planes,  face  normals  a of  electroded  crystal 
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Figure  ia 

Orientation  of  (001)  platelet  for  determining  electric  field 
coefficient  of  0^  through  longitudinal  piezoelectric  effect, 
and  of  02  and  through  transverse  piezoelectric  effect 


The  orientation  of 


together  with  the  associated  Bragg  angles  0^  and  0^. 
these  two  platelets  and  the  location  of  the  two  Bragg  angles  0^  and  0^  is 
illustrated  in  Figures  lb  and  lc.  The  last  column  of  Table  7 gives  the 
corresponding  expressions  for  the  quantity  T obtained  from  the  last  column 
of  Table  2.  If  experimental  values  of  the  electrical  field  coefficient 

* - cot  0^(30^/3E)  are  available  for  all  five  Bragg  angles  0^  (i  = 1,2,... 5) 
described,  the  equations  given  in  the  last  column  of  Table  7 for  M , M2...M,. 
represent  5 equations  for  the  5 unknown  piezoelectric  constants  d^»  d-j2’  ^31* 
d24  and  d15. 

According  to  Table  2 an  additional  independent  measurement  of  the  longi- 
tudinal effect  on  a general  (hk.il)  plane  gives  a linear  combination  of  all  5 
piezoelectric  constants.  Choosing  for  this  purpose  the  (2  2 10)  reflection 
gives  the  entries  listed  on  the  last  line  of  Table  7.  The  platelet  orienta- 
tion and  the  Bragg  angle  0^  are  illustrated  in  Figure  Id.  The  six  expressions 
for  M^,  M2...  listed  in  Table  7 represent  6 equations  for  the  5 unknown 
piezoelectric  constants.  Their  best  values  can  be  obtained  from  a standard 
simultaneous  lfc_.,t  squares  fit. 

For  better  accuracy  additional  platelets  representing  different  values 
of  (0k£) , (h0£)  or  (hkil)  may  be  used  ad  libitum.  It  should  be  noted  that 
for  a given  (hk&)  face  higher  order  and  lower  order  reflections  corresponding 
to  different  Bragg  angles  may  be  measured  and  should  give  the  same  value  for 
the  electric  field  coefficient.  For  example,  the  longitudinal  field  coefficient 
of  all  (OOJl)  reflections  should  be  equal  to  d^  for  any  value  of  Z.  Wnile  the 
lower  order  reflections  002,  004,  006,  008  are  of  lower  intensity  than  the 
0010  reflection,  the  reflection  0 0 20  is  also  relatively  strong  and  may  be 
used  to  obtain  an  additional  independent  value  for  d^. 
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Figure  ib 

Orientation  of  platelet  for  determining  electric  Held 
coefficient  of  0/(  through  trunsveree  piezoelectric  effect 


Electric  field  E 


a 


Figure  Ic 

Orientation  of  platelet  for  determining  electric  field 
coefficient  of  through  transverse  piezoelectric  effect 


Figure  Id 

Orientation  of  (115)  platelet  for  determining  electric  field 
coefficient  of  0t  through  longitudinal  piezoe Lectr I c effect 
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In  order  to  reduce  the  experimental  errors  arising  from  the  specimen 
orientation  it  appears  more  appropriate,  however,  to  use  additional  crystal 
platelets  representing  different  values  of  (OkA) , (h0£) , or  (hkf,) , which 
would  lead  to  different  numerical  coefficients  of  the  piezoelectric  constants 
in  the  equations  for  M^,  M5  and  Mg,  respectively. 

4.  Elastic  and  Thermoelastic  Properties  of  Lead  Potassium  Niobate,  Pb^KNb^O^ 5 
Lead  potassium  niobate,  Pb^KNb^O^^,  occurs  in  the  tungsten  bronze 
structure  and  has  electromechanical  coupling  factors  up  to  73  percent 
(Yamada,  1973).  Since  Yamada  (1973)  also  reported  that  the  temperature 
coefficients  of  the  fundamental  resonance  frequencies  of  crystal  plates  cut 
along  different  directions  have  different  signs,  one  may  expect  the  temperature 
coefficient  to  vanish  for  intermediate  directions.  Provided  that  crystals  of 
sufficient  quality  can  be  prepared,  so  as  to  reduce  the  ultrasonic  attenuation 
to  the  intrinsic  limit  corresponding  to  anharmonic  phonon-phonon  interactions, 
and  provided  this  limit  is  sufficiently  low  so  as  to  be  comparable  to 
ultrasonic  losses  in  a-quartz,  lead  potassium  niobate  would  be  a superior 
substitute  for  a-quartz  in  SAW  devices.  As  reported  above  in  section  2.3,  our 
crystal  growth  efforts  on  this  material  met  with  considerable  difficulties, 
but  we  did  succeed  in  obtaining  several  small  pieces  which  appear  to  be 
predominantly  of  single  domain  type.  In  the  following  we  report  our  prelimi- 
nary data  on  the  elastic  and  thermoelastic  properties  and  compare  the  results 
with  those  recently  published  by  Nakano  and  Yamada  (1975)  and  Yamada  (1975). 

Pb  KNb  CL  c belongs  to  the  orthorhombic  crystal  class  mm2  and  is  ferro- 
2 5 15 

electric,  with  the  spontaneous  polarization  along  the  b axis.  The  Curie 

temperature  is  450°C  (Nakano  and  Yamada,  1975). 

In  Table  8 the  three  lattice  parameters  and  their  associated  linear 
thermal  expansion  coefficients  are  listed.  Based  on  the  present  x-ray  results 
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Table  8.  Lattice  parameters  (A)  and  linear  thermal  expansion  coefficients 

(HfVc)'1)  of  Pb2KNb5015  at  20°C. 


L 

a 

b 

c 

L 

Present 

17.772 

17.972 

7.830 

Yamada* 

17.78 

18.05 

3.917 

I ( 2k) 
l '■ar 

I 

Present 

3.09 

1.06 

3.77 

+ 

Yamada 

1.4 

-0.3 

1.7 

*T.  Yamada,  Appl.  Phys.  Lett.  23,  213  (1973) 
^T.  Yamada,  J.  Appl.  Phys.  46,  2361  (1975) 


and  lattice  constant  data  for  lead  niobate,  our  lattice  parameter  c is 
twice  that  of  Yamada' s (1973)  value.  While  the  lattice  constants  agree 
rather  well,  the  thermal  expansion  coefficients  differ  substantially;  along 
the  b axis  even  the  sign  is  different.  Although  different  experimental 
techniques  were  used  (x-rays  in  the  present  work  and  a dilatometric  technique 
by  Nakano  and  Yamada  (1975))  the  differences  cannot  be  attributed  to  this 
alone,  and  are  more  likely  to  point  to  differences  in  comp'r  .ion,  structure 
or  degree  of  disorder. 

The  calculated  and  measured  density  values  in  Table  9 show  rather  small 
differences  between  Yamada' s and  the  present  results,  perhaps  too  small  to 
point  toward  substantial  compositional  or  structural  differences.  However, 
in  order  to  explore  the  effect  of  ordering  phenomena  on  property  differences, 
a considerable  amount  of  further  vnrk  on  possible  correlations  between 


3 

Table  9.  Density  values  (in  g/cm  ) of  Pb2KNb^01^ 


Calculated  (X-Ray) 

Measured 

Present 

6.15 

6.23 

Yamada* 

6.12 

6.14 

Yaraada,  Appl.  Phys.  Lett.  22_,  213  (1973) 


thermal  history  of  the  samples,  their  relative  x-ray  intensities  and  their 
elastic  and  thermoelastic  properties  is  required.  Since  the  measured 
densities  are  larger  than  the  calculated  values,  the  crystals  used  both  by 
Yamada  and  in  the  present  work  may  have  had  a small  excess  of  lead. 

In  Table  10  the  temperature  coefficients  of  the  transit  times  X for 
longitudinal  waves  along  the  three  crystallographic  axes  are  compared  with 
Yamada' s (1975)  values  for  the  temperature  coefficients  of  the  corresponding 
principal  resonance  frequencies  f.  The  relation  between  these  two  quantities 
is  (Chang  and  Barsch,  1976) 


1 1 

[3f] 

_ _ i 

2kZ 

f 1 

T 

2 2 2 
Z -k  (1-k  ) ^ 

Here  k denotes  the  electromechanical  coupling  factor,  and  L the  solution  of 
the  equation 

V 

tan  Z = Z/k^ 

For  small  values  of  k and/or  (9k/9T)  the  temperature  coefficients  of  f and 
T should  be  numerically  equal,  but  opposite  in  sign.  According  to  Yamada 
(1975)  the  X-  and  Y-cut  specimens  have  zero  or  smaxl  values  of  k for  longi- 
tudinal waves,  respectively,  so  that,  unless  (9k/9T)  is  excessively  large, 
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Table  10.  Temperature  coefficients  of  transit  time  for 
longitudinal  waves,  (l/l) (9t/9T) , and  of 
fundamental  resonance  frequency,  (1/f) (3f /3T)  , 
(.both  in  10~6(°C)-1)  for  X,  Y and  Z cut  of 
PboKNb,0.  _ at  25°C. 

I J 1 J 


X-cut 
(J_  c-axis) 

Y-cut 
(_[_  a-axis) 

Z-cut 
(J_  b-axis) 

(1/x) (9t/9T) 
(Present  Work) 

-3.03 

- - 

84-00 

74.82 

(l/f)(9f/3T) 

(Yamada*) 

L 

200 

25 

-30 

*T.  Yamada,  J.  Appl.  Phys.  _46,  2894  (1975) 


(1/t)(9t/9T)  ~-(l/f)Of/9T).  The  data  in  Table  10  indicate  that  this 
relation  does  not  at  all  hold.  Again,  this  may  be  attributed  to  differences 
in  composition,  structure  or  degree  of  order. 

The  orientation  of  the  cartesian  x,  y and  z-coordinate  axes  along  the 
crystallographic  c,  a and  b-axes,  respectively,  follows  the  convention  of 
Yamada  (1975).  These  directions  were  determined  optically  by  means  of  the 
quartz  wedge  method  (see,  e.g.  Hartshorne  and  Stuart,  1964).  In  this  manner 
a value  of  An  = n&  - n^  = 0.064  for  the  difference  of  the  refractive  indices 
along  the  a and  b axes  was  determined  at  3G0"'K,  as  compared  with  a value  of 
An  - 0.055  given  by  Nakano  and  Yamada  (1975).  The  a and  b axes  can  be  dis- 
tinguished on  the  basis  of  the  sign  of  An  and  the  fact  that  according  to 
Nakano  and  Yamada  (1976)  the  a-axis  is  the  slow,  and  the  b-axis  the  fast 


direction. 
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In  Table  11  the  on-diagonal  elastic  constants  and  their  temperature 
coefficients  are  listed.  The  assignment  of  the  coordinate  axes  follows 
Yamada  (1975),  that  is  x||c,  Y | |a,  and  z||b.  Except  for  c55  considerable 
differences  in  the  elastic  constants  of  Yamada  and  the  present  data  occur. 

No  other  explanation  except  uhe  above  mentioned  conjecture  about  differences 
in  the  nature  of  the  sample  can  be  offered. 

It  should  be  noted  here  that  the  as-grown  boules  obtained  by  us  are  of 
considerably  higher  quality  than  those  of  Nakano  and  Yamada  (1975).  Whereas 
the  specimen  on  which  the  above  elastic  data  were  obtained  do  net  show  a 
significant  amount  of  elastic  twinning  under  a polarizing  microscope  (i.e. 

90°  domains  for  the  spontaneous  polarization  Pg||b-axis,  or  presence  of  domains 
with  their  a and  b axes  interchanged)  the  specimen  obtained  by  Nakano  and 
Yamada  (1975)  were  thoroughly  twinned  and  had  to  be  de-twinned  by  strain 
annealing  in  an  electric  field.  It  is  conceivable  that  by  this  procedure 
residual  strains  are  introduced  which  could  affect  the  elastic  constant  data. 

Although  the  specimen  used  in  the  present  investigations  were  virtually 
free  of  elastic  twins,  the  presence  of  electrical  twins  (i.e.  180°  domains 
for  the  spontaneous  polarization)  cannot  be  ruled  out,  since  this  type  of  twin 
cannot  be  detected  under  a polarizing  microscope.  In  such  twins  all  five 
piezoelectric  constants  of  the  class  mm2  (that  is,  e15’  e24’  e31’  °32’  °33) 
show  opposit0  signs.  However,  since  the  piezoelectric  stiffening  term  in  the 
Kristof f el  tensor  depends  only  on  the  square  of  the  piezoelectric  constants 
the  results  for  the  elastic  constants  are  not  affected  by  the  presence  of 
180°  domains. 

According  to  Table  11  the  temperature  coefficient  of  c^  is  positive,  so  that 
one  may  expect  temperature  compensated  cuts  to  exist  for  longitudinal  waves  along 
a direction  close  to  the  X-direction  (c-axis) . Unfortunately,  this  is  not  where  the 


D 12  2 

Table  11.  Adiabatic  on-diagonal  elastic  constants  c (10  dynes/cm  ), 
and  their  isobaric  temperature  derivatives  (9cjV/  T) 

pV  p 


(107  dynes/cm^)  for  Pb^KNb^O^^  at  25  C 


yv 

11 

22 

33 

44 

55* 

66 

- - 

D 

c 

yv 

Present 

1.46 

1.03 

1.45 

.450 

.600 

Yamada 

1.66 

1.63 

1.91 

0.66 

0.57 

0.63 

r D 

c 

yv 

Present 

— 
+ .465 

-18.39 

-29.49 

-2.83 

-5.50 

L 9T  J 

P 

This  constant  could  not  be  measured  since  with  the  presently  available 
crystal  no  echoes  of  the  ultrasonic  signal  could  be  received. 

^T.  Yamada,  J.  Appl.  Phys.  4^6,  2894  (1975) 


electromechanical  coupling  factor  reaches  its  maximum  (Yamada,  1975).  However 
for  the  exact  determination  of  the  orientation  of  temperature  compensated  cuts 
for  bulk  or  surface  waves  accurate  values  of  all  elastic,  dielectric  and 
piezoelectric  constants  and  their  temperature  derivatives  are  needed.  It  is 
conceivable  that  the  electromechanical  coupling  factor  for  these  directions 
turns  out  to  be  larger  than  for  a -quartz  and  a-berlinite  (Chang  and  Barsch, 
1976). 

Our  crystal  growth  efforts  to  obtain  still  larger  single  crystals 
required  for  the  measurement  of  the  remaining  elastic  and  thermoelastic 
constants,  and  for  the  independent  measurement  of  the  dielectric  and 
piezoelectric  constants  are  being  continued. 
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